NOTES ON THE DRUDE MODEL 



KYLE MCELROY 



Abstract. Here are some notes about the classical theory of metals and the 
Drude model. The classical theory of electrons in metals is used to describe 
the basics of metallic behavior. These include their electrical conductivity, 
heat conductivity, reflection of optical wavelengths, etc. 



1. Assumptions and basics 

The Drude model was developed at the turn of the 20th century by Paul Drude. 
It came a few years after J.J. Thompson discovered the electron in 1897. It predates 
quantum theory, but still can tell us a lot about electrons in metals. 

As background the the model we should get to know the electrons and how many 
we are dealing with. We're going to keep the valence assumption. This assumption 
rests on the intuition that the core electrons will be more tightly bound to their 
nuclei and hence will not be free to wander around and contribute to conduction. 
Essentially this lowers the number of electrons from Z to Z c where Z c is the number 
of conduction electrons. 

So in a sample of metal say Sodium (Na) the density of conduction electrons, n 
is: 
(1.1) 

n = Na ^ = 6 .02 x 10^ato mS /mol le/at ° m -\ XW69/m3 = 2 x 10 28 e/m 3 
A ' 29g/mol ' 

where Na is Avogadro's number, p m is the density of the metal, A is the atomic 
number of the element and the numbers are for Na. 

For the actual model we're going to eliminate all the electron ion interactions 
and replace them by a single parameter. 

(1) We will treat collisions between e's and ions are instantaneous, uncorrclatcd 
events. 

(2) We will ignore all other interactions (i.e. potentials from ions or other 
electrons) except applied fields. This means that electrons travel in straight 
lines between scattering events. 

(3) Probability of an electron having a collision in a time interval dt is dt/r. 
And t does not depend on the electron position or momentum. 

(4) Collisions 'thcrmalize' electrons. This means that after a collision the elec- 
trons have the temperature of the local environment. 

2. Equations of motion 

The first thing you need is to figure out how an electron's momentum, on average, 
will evolve over time. To do this we'll just find the average equation of motion for 
an electron. To find this let's start with the momentum of an electron at time t, 
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p(t), and find it at time t + dt. If the electrons had a collision it would on average 
have no momentum (j>c(t + dt) = 0) at time t + dt and by the third assumption 
above this has the probability, P c = dtjr . This means that the probability of no 
collision is P nc = (1 — dt/r) this is because P c + P nc — 1. If there were no collision 
the electrons would have evolved normally and the electrons momentum becomes, 
Pnc(t + dt) = p(t) + F(t)dt. This makes the new momentum: 

(2.1) p(t + dt) = P c ■ p c (t + dt) + P nc ■ p nc (t + dt) = (1 )[p(t) + F(t)dt] 

T 

Using this to find the derivative take 

(2 2) m = p\t + dt)-p\t) = p(t) p 

[ ' ' dt dt t U 

And you have the equation of motion (EoM) averaged over electrons. 
Of note there arc a few regimes and solutions to consider. 

• If F(t) = the solution to this homogeneous equation is p(t) = p(0)e _ */ T 
which is why r is called the relaxation time. If you impart momentum to 
the electrons on average they will relax back to no momentum exponentially 
with a time constant r 

• With a constant F you can show that the solution to the momentum p(t) = 
p(0)e-*/ T + Ft 

• After a long time, t>T, the exponential term becomes negligible leaving 
p(t) = Ft 

3. Ohm's law 
You all know Ohm's law: 

(3.1) V = IR 

with V is the Voltage applied to a metal, / is the resulting current, and R is the 
proportionality constant. The main empirical fact here being the proportionality 
of current to the applied Voltage. This is what we're going to try to predict. But 
let's recast it in a form that is not dependent on the geometry of the experiment: 

(3.2) J=aE 

where a is the conductivity of the metal and j is the current density. Since j is 
the current density it is the number of electrons passing a given point or j = l/nev 
where n is still the electron density, e is still the electron charge, and v is the average 
drift velocity of the electrons. 

In an applied electric field E the EoM for long times, t ^> r, gives us: 

p(t) = eEr 

_. . eET 

v{t) = 

m 

. Plugging this into the expression for the current density gives us 

(3.3) J=^1S 

m 

which is Ohm's law with a = ne 2 T /m. 

There are a couple of implications of this. For a metal like Na with a resistivity, 
PNa = 1/c = 50 nil-m the relaxation time is about 10~ 14 s. 



4. AC CONDUCTIVITY AND THE SHINY 

In HW^l you proved that a complex form of Ohm's law still holds for a metal in 
the Drude model with an oscillating electric field incident on a metal (i.e. E(t) = 
Re{E(aj)e tu ' t }). With an E-ficld of this form Ohm's law takes the form 

-», . ne 2 T 1 -±, , 

M = — - — -r^H 

m illjt + 1 
where the u> dependent conductivity is given by 

ne 2 T 1 
o-(w) 



m iujT + 1 



Since the electric field is oscillating electric we have an electromagnetic wave. So 
we should look at Maxwell's equations which with no net charge look like: 

V -E = V -8 = 

dB - - , ff 

Vx£=-- V x B = noj + Moeo^f 

Since we have the complex Ohm's law we can replace the current density in the 
last equation and simplify it to be 

V x B = (yUocr(w) + iuj[j,oeo)E 

Let's take the curl of the third equation and use the first one and what we just 
found to simplify it to 

-> 9 •* dV x B 
VxVx£= —vE = — = iuj{no<j{tjj) + icu^o^ojE 

After some rearranging this can be written as 

_S9 ^ w 2 ,i<j(u)) . j» 
-V 2 E= — (— y —^ -\)E 
c z wen 

Now notice that if you take the limit where lot 3> 1, <r{u>) — > gqJilot leaving us 
with 

2 2 

-V 2 E=^( 2 -l£ 

Now if you assume that the electric field has a plane wave spatial structure (i.e. 
E(t) = Re{E(cu)e luJt e ik ' r }}) then the above equation simplifies to the dispersion 
relation 

k 2 c 2 = (uj 2 - UJ 2 p ) 

with uip = ne 2 /meo being the plasma frequency. 

This dispersion relationship tells us that k is imaginary for uj < lu p which is the 
result for an exponentially decaying electric field. This means that the EM wave 
doesn't propagate into a metal at these lower frequencies and instead is reflected. 
On the other hand for u> > lo p this gives a real k which is the form for a travelling 
wave. This implies that for frequencies higher than the plasma frequency metals 
become transparent. Which experiments verify (see figure ). 

For metals the plasma frequency is in the ultra violet and for Na it is about 
uip = 10 15 s _1 which corresponds to a wavelength of light at A 200nm. 
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Note: We disregard the direct magnetic field effects because they are much 
smaller than the electric field forces. 

5. Joule heating 
See the solution set for HW#1 

6. Hall effect 

We worked out what happened with a constant electric field put onto a metal. 
But what if you use a Voltage to drive current through a metal and put a perpen- 
dicular magnet field on it. 

7. Thermal conductivity 

The electrons in a metal aren't just good conductors of electricity they are also 
good at conducting heat. So we're asking about the amount of heat that travels 
along a sample of metal as heat is applied. There is a simple phcnomcnological 
model that is analogous to Ohm's law, it is called Fourier's law and it looks like: 

f = -kVT 

The minus sign is because heat flows to lower temperature, and j q is the heat 
current density AQ/Time-A, heat per time per area in units of [Watts/m 2 ] and n is 
the thermal conductivity which has units of [Watts/mK] . So what does the Drudc 
model predict? 

A couple of assumptions are important here. First we need to remember the 
last assumption of the Drude model. This assumption means that after a collision 
an electron carries the thermal energy of the local environment. Also we need to 
assume that T varies a little over I (the mean free path electrons travel before 
scattering) . 

To calculate the heat delivered to a point we need to simplify this to a 1-D 
problem and then we need to calculate the heat it gets from the left and subtract 
the heat it gets from the right. Let's write the thermal energy an electron at 
temperature T has as £(T). Then at a point x the average electron coming from the 
left brings with it an energy £(T[x — vt]) the average electron from the right delivers 
£(T[x + vt]). So for electrons of density n and average velocity v (remembering 
that half will travel towards the point x and half away) this leaves us with 

j q = -nv{£(T[x — vt] — £(T[x + vt]} 

1 . d£ . , dT . „ 

al ax 
Now note a few things n4§* — yM< — hipf = c v or the specific heat per volume. 
Also when we generalize to 3-D we note that < v 2 . >=< v 2 >~<v 2 z >= | <w> 2 . 
Putting all this together we have 

So the prefactor up there is equal to the thermal conductivity. 



There is another way to look at the specific heat. It's called the Wiedemann- 
Franz ratio. Which is given by 



k m < v > 2 TCy 



oT 3e 2 

Now in a classical gas using the equipartition theorem the total thermal energy per 
volume is given by E = n^kgT so the specific heat per volume (the derivative of 
the energy with respect to energy) is cy — n^ks- And also note that the total 
thermal energy is the kinetic energy is given by \m < v > 2 — n^ksT which we can 
use to get rid of the in < v > 2 in the W-F ratio. All together this gives 

^ = |(-) 2 = l-HxlO- 8 (^/X 2 ) 

This number only has fundamental constants in it and it is 1/2 the real value you 
get experimentally. But is really close for such a simple model. 



